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Abstract 



We study the universal properties of distributions of eigenvalues of random matrices in 
the large N limit. The distributions fall in universality classes characterized entirely by 
the support of the spectral density. 
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1 The theorems 



Random matrices have found a wide range of applications in solid state physics, nuclear 
physics and high energy physics. An important observable in most of these applications is 
the correlator between eigenvalues of the random matrices. In addition to being random, 
the matrices are often large, and it is natural to take the so called large N limit in many 
of the applications. The general set up is as follows: we consider an ensemble of random 
matrices 

P (0) = i e -^^), Z = /# e" mv W, V{<j>) = ]T f ^ (1) 

z j k 

and in this ensemble we can ask for expectation values of a certain observables : 

</(0)> = |#P(^)/(0). (2) 

In the following we will assume that the random matrices <f> are N xN Hermitian matrices, 
but similar results can be proven for ensembles of complex matrices matrices. 

The one- and two-point "resolvents" are defined by 

G{z u z 2 ) = / Tr ^— Tr^-\-/Tr^-\/Tr^— \, (4) 

\ z\-(p z 2 -(f>/ \ zi — <p/\ z 2 -(p/ 

They are related to the spectral density and the correlator of densities 

1 N 

p(X) = {S(X)), 5(A) = -X;*(A-A<) (5) 

i=i 

and 

p(X,X') = (S(X)S(X'))-(S(X))(S(X')) (6) 

in the following way 

P(A) = 2^( G(A + i£) " G(A "^ ) ) (7) 

and 

p(A, A') = ^ (^) 2 (G(+, +) + G(-, -) - G(+, -) - G(- +)) , (8) 
where G(±, ±) = G(X ± ie, X' ± ie')- 

The following theorem was proven in |jj (and partly rediscovered by the authors of [Q], 
after whom it is sometimes called Brezin-Zee universality): 

Theorem 1: Assume that the spectral density p(X) has support in a single interval 
in the large N limit. Then the two-point correlator (and thus p(A, A')) is universal, 
i.e. independent of V((f>), and given by (in the case of a symmetric potential, such 
that the support is [—a, a]) 

G(zuz 2 ) = Tl 1 L 2+ (^ 2 ) + fe 2 -a 2 



(9) 
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For the general formula for an arbitrary potential, as well as for complex matrices, rather 
than Hermitian matrices, we refer to O. In fact much more was proven in jjj]. Any 
multi-point resolvent 

G{z u ... lZn ) = iV 2 "- 2 (-j-Tr ■■■^Tr - J— \ , (10) 

\iV iV z n -0/ con „ 

where conn refers to the connected part of the multi-point resolvent, is also universal. 
Closed and very simple expressions were given for the multi-point resolvents in the case 
where matrix ensemble consisted of complex matrices. Note that the generalization of (||) 
is 

(S(X 1 )...S(X n )) conn = j^(^) n £ (-l)^G{z l +a 1 ie,...,z n + a n ie). (11) 

ffc=±l 



In addition it is possible to develop a systematic 1/N 2 expansion: 

1 



oo 



G(z 1 ,...,z n ) = ^2 Tf2h G h (zi,...,z„,), (12) 



h=0 
oo -. 

F = HjpK F h- (13) 

h=0 

Again the two-point resolvents and in this case even the corrections to the spectral density 
p(X) as well as the corrections to the free energy N 2 F = \ogZ are all universal |J. These 
considerations have later been extended to supermatrices 

In the above discussion universality means that the resolvents can be presented in a form 
independent of the potential. Of course the function (H) depends on the potential through 
the endpoints of the eigenvalue distribution p(\). The same remark is true for the multi- 
point resolvents and for higher order expansions in 1/N 2 . However, for the n-point resol- 
vent calculated to non-trivial order h in the 1/N 2 expansion only 2(3/i — 2 + n) additional 
parameters \ are involved for an arbitrary potential (see def. (|27|)). Yet, some of these 
might vanish, as will be the case for a Gaussian potential where all additional parameters 
vanish but one. The somewhat misleading statement found in Q that all multi-point 
density correlators vanish to order 1/N n is a trivial fact which follows from the factoriza- 
tion property of the large N expansion. The first non-trivial term is of order l/N 2n ~ 2 for 



the n-point density correlator (11) and it is explicitly non-Gaussian in nature for n > 2 
because of the additional parameters. These observations can be summarized as follows 

Theorem 2: If the spectral density in the large N limit has support [—a, a] all higher 
1/A^ 2 corrections to the free energy F, to the n-point resolvents, to the spectral 
density and any multiple eigenvalue correlators are universal. The correction to the 
n-points resolvents are of the form 

R h ({z 2 -a 2 },{MP},a) 
G h ( Zl ,...,z n ) = V , J , (14) 

lm=i(z 2 -a 2 ) 



where the Rh are simple rational functions of the arguments, which can be deter- 
mined iteratively in the 1/A^ 2 expansion. 
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We refer to Q for details, generalizations to the situation where the support is not sym- 
metric and to || for complex matrix models. 

All these considerations have been based on the assumption that p(\) has support in one 
interval on the real axis. More general situations can occur for complicated potentials, 
and in fact it is natural for applications in solid state physics to have such a situation 
since the eigenvalues of <p are viewed as (part of) the eigenvalues of a Hamiltonian which 
can have a band structure. In this paper we prove 

Theorem 3: Assume that the support of the spectral density p(\) in the large 
N limit consists of s intervals \x2j1 %2j-i]i j = 1, • • • ,s. The large N limit of the 
two-point resolvents Go (21,22) (and therefore the correlators of eigenvalues) fall in 
universality classes characterized by s and they depend only on the potential V 
through the endpoints X{. For each s > 1 it is possible to generalize the statements 
in Theorem 1 and 2 concerning universality. 



Let us here present the explicit solution, analogous to @), in the case where the support 
of p(A) consists of two intervals [3:4,3:3] and [x2,xi], X4 < X3 < X2 < x%, and where we for 
simplicity assume that V(4>) is an even function of <f> such that X4 = —x\ and X3 = — x 2 . 
If we denote x±, X2 by a, b we have 

1 / o , (2f-a 2 )(2l-6 2 ) + (2 1 2 -6 2 )(2l 



G (2i,2 2 ) = — so -2 + 



4(21 - 2 2 ) 2 ^ ' ^ z 2 _ a 2^ z 2 _ b 2^ z 2 _ ^2 _ 6 2) 

+ {a + bfm 1 

4 K ( k ) y/(zl - a 2 )(zf - tf){z 2 2 - a 2 )(4 - b 2 ) ' 

where k = 2y/ab/(a + b), and where E(k) and K(k) denote the complete elliptic integrals 
of first and second kind. 

In the next section we will briefly outline how to prove theorem 3. 

2 Proofs 

The basic tool for proving the above statements is the loop equation in the following form 
(1): 

JCG(z) = G 2 (z) + 1 ^G(z,z). (16) 
In this formula fC is a linear operator 

(*/)« - / a?) 

Jc 2m z — to 

and the contour encloses all singularities of G(z), but not z. Note that the support cr(A) 
of p(X) coincides with the singularities of G(z). In large limit they will be located at 
s cuts on the real axis. If the Xj, i = 1, ... ,2s denote the positions of the endpoints of 
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the cuts [x\ > X2 > • • •) the solution to ([16]) in the large N limit can be written in closed 
form since we can discard the term involving G(z, z). We get 



G (z) 



1 r du V'{u) 



2 Jr 2ni. : - w\l JJ 



2s 



Z Xj, 



UJ — Xi 



(18) 



In addition the endpoints of the cuts are uniquely determined by the fact that Gq(z) ~ l/z 
for large \z\, and a certain stability requirement ||. The first requirement leads to s + 1 
equations: 

i r dio v'{uj)uj k 



2j c 2m ^ n 2. i(w _ x . 
If we define the polynomial M{z) by 

M(z) 



Sk,. 



k = 0, . . . , s. 



(19) 



(20) 



where Coo is a contour at infinity in the complex plane, the requirement of stability implies 
@ ' 



d\M(\) A 



X2k+1 



2s 



\ i=l 



0, fc = l,...,s-l. 



(21) 



It is this kind of boundary condition which leads to the appearance of hyperelliptic integrals 
in formulas like (ITf 



Note that it follows from (18) that the planar limit of the spectral density is 

,) , Aecr(A). 
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p(A) = — |M(A)| 



\ 



2* 



i=l 



(22) 



Having obtained the complete solution Gq(z) in the large N limit one can solve the loop 
equation (16) iteratively as an expansion in 1/iV 2 , and at a given order h it is possible to 



construct the multi-point resolvent Gh(z\, . . . , z n ) from Gh{z). The key ingredient in this 
construction is the so called loop insertion operator 



d 



dV{z) 



-E 



k d 



yk + l 



dgk 



(23) 



By definition it follows thatj^j 
G(z 1 ,...,z n ) 



d 



d 



In particular, we have 



dV{z n ) dV{z x ) 



G(zi,z 2 ) 



dV(z 2 ) 



F, N 2 F = log Z, n>2. 



G{z x ) 



(24) 



(25) 



2 For G(z) l/z has to be added on the r.h.s. 
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which allows us to construct Gq(zi, Z2) from (|1S|). The steps are in principle elementary, 
but a considerable amount of algebra is needed in order to prove that Gq(z\, Z2) is universal, 
and of the form (15) for the example of two cuts with an even potential V(z). The details 
will be published elsewhere Q. Let us only mention here that the main complication 
compared to the single cut case (i.e. s = 1) is that the kernel of fC — 2Gq{z) is not zero 
but given by 

z l 

Ker(/C - 2G (z)) = Spanj I = 0, . . . , , s - 2}, (26) 

after taking into account the asymptotic of G(z). 

The s-cut solution will be characterized by 2s classes of so called moments, given by 

Following the derivation for s = 1, but with considerable algebraic complications, one then 
can show by iteration that Gh{z\, . . . , z n ) can be written in the form 



. . . ,Z n ) — , (2»J 

/nfiilT™=i(^-^) 



where i?^ is a rational function of the poles and the moments with 1 < I < 3h — 2 + n. 
However, contrary to the case s = 1 the coefficients in are no longer rational functions 
of the Xi, but involve hyperelliptic integrals as well. 



3 Comments 

The universal features of the single cut solution of the random matrix models generalize 
to an arbitrary number of cuts. As was shown in ||, 10] the universality of the one-cut 



solution is deeply connected to the theory of intersection indices on the moduli space of 
punctured Riemann surfaces. It is unknown if a similar mathematical underlying structure 
can be related to the multi-cut solution. The fact the the coefficients in (|28l ) can be 
expressed in terms of hyperelliptic integrals gives some hope in this direction, since similar 
hyperelliptic integrals are encountered for the random matrix model coupled to 0(n) spins 
and this model is known to be related to integrable hierarchies. 



References 

[1] J. Ambj0rn, J. Jurkiewicz and Y. Makeenko, Phys.Lett B251 (1990) 517. 
[2] E. Brezin and A. Zee, Nucl.Phys. B402 (1993) 613. 

[3] J. Ambj0rn, L. Chekhov, C.F. Kristjansen and Yu. Makeenko, Nucl.Phys. B404 (1993) 
127. 

[4] J.C. Plefka, Nucl.Phys. B444 (1995) 333; B448 (1995) 355. 

[5] J. Ambj0rn, C.F. Kristjansen and Yu. Makeenko, Mod.Phys.Lett. A7 (1992) 3187. 



6 



[6] F. David , Nucl.Phys. B348, (1991) 507. 
[7] J. Jurkiewicz, Phys.Lett. B245 (1990) 178. 

[8] G. Akemann, Higher genus correlators for the hermitian matrix model with multiple 
cuts, NBI-HE-96-30, IPT-UH-08/96, frep-th/9606004 . 

[9] L. Chekhov, Geometry and Physics, 12 (1993) 153. 

[10] L. Chekhov, Matrix models and geometry of moduli spaces, Steklov Mathematical 



Institute Preprint, [hep-th/9509"001 



7 



